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NUCLEAR  FUSION  LN  A  SOLID  VIA  A  BOSE  BLOCH  CONDENSATE 
Introduction 

This  paper  discusses  a  concept  of  solid  state  fusion  based 
on  the  formation  of  a  D’  boson  Bloch  condensate  (BBC)1  within  a 
palladium  deuteride  host  lattice.  Evidence  for  such  solid  state 
fusion  has  been  presented  by  Fleischmann  and  Pons2,  who  have 
described  episodic  generation  of  anomalous  heat  in  Pd  cathodes 
following  extended  overvoltage  electrolysis  of  LiOD  solutions. 
Other  evidence  for  solid  state  nuclear  interactions  in  PdDx 
includes  the  observations  of  surface  hydrogen  isotopic  anomalies 
in  Pd  cathodes  subsequent  to  such  electrolysis3-4  and  the  possible 
observation  of  a  small  number  of  20+  Mev  energy  release  events  in 
a  Si  charged-particle  detector  adjacent  to  a  Pd  electrode  charged 
with  D'  ions  by  ion  implantation5. 

The  bosonic  Bloch  condensate 

The  D’  BBC  is  a  collective  bosonic  condensate  derived  from 
weakly  bound  D’  ions  which  are  well-screened  but  constrained  to 
an  approximately  periodic  arrangement  through  electrostatic 
interaction  with  a  host  lattice.  The  resulting  Bloch  symmetry 
permits  the  condensate  to  exist  at  low  boson  concentration  c  = 
Ng/ft,  where  N,  is  the  number  of  condensate  bosons  and  NL  is  the 
number  of  unit  cells  in  the  host  crystal.  The  BBC  is  a 
collective,  cooperative  entity  which  contribut  fractional 
charge  per  unit  cell.  It  undergoes  integt  occupation 
fluctuations  at  individual  sites,  expressed  by  Wannier  functions. 

Palladium  hydride  and  conditions  favoring  D*  BBC  formation 

The  PdH„  and  PdD„  systems  have  been  studied  extensively6. 
PdD,  differs  from  PdH„  by  having  a  lower  equilibrium  pressure  at 
any  given  x.  Both  materials  are  superconducting  at  low 
temperature  with  PdD„  having  a  higher  transition  temperature. 
This  higher  transition  temperature,  and  the  higher  diffusivity  of 
D  as  opposed  to  H  in  Pd  are  contrary  to  expectations  based  on 
increased  mass.  Absorption  of  D  into  PdD,  expands  the  lattice 

and  causes  locai  strain  when  x  is  smell.  This  strain  energy  is 

Manuscript  approved  January  12.  1990 


1 


a  component  of  the  chemical  potential,  which  is  very  positive  at 
low  x,  becoming  very  negative  as  x-»l.  As  x-*l,  D  additions  remove 
local  strain  due  to  isolated  vacancies.  This  x-behavior  of  the 
chemical  potential  suggests  that  values  of  x  very  near  x=l  should 
be  favorable  for  formation  of  a  D+  BBC.  Very  near  x~l,  the 
limited  availability  of  vacancies  dictates  that  some  unit  cells 
will  contain  more  that  1  D+  ion.  Without  formation  of  a  0*  BBC, 
these  sites  must  be  doubly  occupied,  creating  high  local  lattice 
strain,  hence  a  high  positive  chemical  potential.  When  a  D'  BBC 
forms,  the  excess  D+  charge  is  shared  by  all  periodically 
equivalent  sites  in  the  crystal,  avoiding  local  lattice  strain. 
Hence  the  D+  BBC  should  be  thermodynamically  favored. 

D+  BBC  fusion  reactions 

Nuclear  reactions  from  a  D+  BBC  state  in  a  crystal  lattice 
are  possible  because,  from  such  an  initial  state,  the 
electrostatic  barrier  problem  of  free  space  nuclear  interaction 
may  be  eliminated.  Fusion  is  a  result  of  three  D*  BBC  properties: 
screening  by  the  lattice  for  c  <<  1,  the  resulting  large  zero 
point  motions  associated  with  this  screening,  and  the  multiple 
boson  occupation  of  single  lattice  sites.  {Here  BBC-D*  refers  to 
the  D+  ions  that  have  become  part  of  the  BBC,  and  does  not  include 
D+  ions  which  are  part  of  the  host  lattice.  Henceforth  D*  shall 
refer  to  BBC-D+  unless  otherwise  indicated.)  The  large  zero 
point  motion  spreads  out  the  D+  charge  density  within  a  unit  cell  . 
The  resultant  reduction  in  the  self-interaction  associated  with 
multiple  occupation  virtual  states  increases  the  time  constant 
for  electrostatic  interaction  relative  to  that  for  nuclear 
interaction,  leading  to  a  separability  of  interaction  modes,  as 
discussed  further  below.  Since  each  D+  BBC  ion  "sees"  an 
approximately  periodic  potential,  the  single  particle  density  at 
a  given  site  is  derived  from  a  superposition  of  different  Bloch 
symmetric  eigenstates.  The  associated  density  at  a  given  site 
can  be  viewed  as  either  a  superposition  of  non-stationary , 
Wannier  state  densities,  each  of  which  is  localized  primarily 
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around  an  individual  site,  or  as  a  superposition  of  time 
independent  Bloch  state  densities,  each  of  which  extends 
throughout  the  lattice.  In  the  evaluation  of  the  transition 
matrix  element  associated  with  fusion,  the  initial  state  wave 
functions  and  fusion- inducing  changes  in  the  potential  result 
from  non-number  conserving,  many-body  fluctuations  associated 
with  multiple  occupation  of  Wannier  states  at  a  given  site.  Site 
multiple  occupation  causes  deuteron  wave  functions  to  overlap, 
inducing  fusion. 

BBC  formation  using  bosons  which  interact  only  with  the  lattice 

We  first  consider  a  BBC  formed  from  autually-non-interacting 
bosons.  These  bosons  interact  only  with  a  host  metal  lattice  and 
its  itinerant  electron  population.  Boson-boson  interactions  are 
ignored . 

We  consider  the  case  of  the  H,  bosons  distributed  within  a 
crystal  containing  host  unit  cells.  It  is  assumed  that  the 
bosons  reside  in  potential  wells  within  the  unit  cell  of  the  host 
metal.  Fractional  charge  is  allowed  at  each  site  so  that  the 
potential  provided  by  the  lattice  is  periodic  over  a  large  enough 
scale  that  Bloch  functions  apply.  The  single  particle  wave 
functions  ^Bto<h  (k,  r)  for  the  ideal  lattice  have  the  property 

-  sl;o<h(k,r)  exp(ik.R)  .  (1) 

The  fluctuation  properties  are  shown  by  expanding  in  terms  of 
Wannier  functions 

Nl 

^Bloeh  (  fc*  r  )  exp(-«  (k)  t/ri)  «  (1/Nl)1/2  E  *.(r,t)  exp(ik*RJ  ,  (2) 

s=l 

where  R  and  R,  are  Bravais  lattice  vectors,  «  (k)  is  the  band 
energy  of  <^Bloch,  an<*  *  is  the  crystal  momentum.  Each  is  a 
Wannier  state,  which,  in  the  harmonic  approximation,  is  suitably 
approximated  at  t=0  by  the  ground  state  wave  function  of  a 
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parabolic  well, 

<Mr)  =  (2/ wa2) 3/*  exp(-  rs2/a2)  ,  (3) 

where  rs  =  r-R*,  and  a  is  the  classical  turning  point  of  the  well. 

The  many  body  wave  function  associated  with  the  occupation 
of  Nb  band  states  (denoted  by  eigenvalues  «p)  possessing  NB 
coordinate  dependencies  r  is  given  by 

Mb_ 

V(«p,r)  =  (l/N  B!  )  1/2  S  ||*  Bloch  (  /  rm)  /  «> 

{ rm}  m=l 


where  the  sum  over  {rm}  includes  interchange  of  each  coordinate 
rm  with  the  remaining  NB-1  coordinates,  ensuring  that  *  is 
suitably  Bose  symmetric.  Substituting  Eq.  2  into  Eq.  4,  we  find 
that 

N  /  2  ^ 

¥(cp,r)  =  (  1/Nb  ! ) 1/2  Z  (1/Nl)B  {  M  Z%.(r0)  expUVRJ  >  •  (5) 

{rm}  m=ls=l 

In  the  bracketed  product  of  summed  terms  in  Eq.  5  there 
exist  subsets  of  terms  containing  multiple  values  of  m 
contributing  to  fixed  site  s  Wannier  functions  <^,(rB)  .  All  terms 
for  which  n  values  of  m  contribute  to  any  of  the  ^s(rm)  correspond 
to  n-fold  occupation  of  the  lattice  site  s.  For  large  NB,  small 
c,  and  small  occupation  number  n  there  exist  (NL)ra  cn/nl  terms 
corresponding  to  n-fold  occupations  by  distinguishable  bosons. 
For  indistinguishable  bosons,  these  terms  become  n! -degenerate 
due  to  the  equivalence  of  n!  permutations  of  the  coordinate 
dependencies  in  each  of  these  terms.  Thus  there  are  ( NL ) ^  cn  terms 
corresponding  to  n-fold  occupations  by  indistinguishable  bosons. 
For  notational  purposes,  we  regroup  terms  in  Eq.  5  that 
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correspond  to  different  n-fold  occupations.  The  corresponding 
collection  of  terms  we  designate  as  #(n,k^,  rj.  Thus 

N. 

*  ( « j,.  r)  -  £  *  (n.  r) 

n=l 

In  the  next  section  we  will  he  concerned  with  4-fold  occupations, 
which  we  call  quad-bosons. 

Nuclear  reactions 

We  now  consider  the  properties  of  O'  BBC  bosons  and  the 
resulting  nuclear  interactions  that  are  implied  by  self¬ 
interaction.  As  discussed  in  the  next  section,  the  magnitude  of 
the  electrostatic  sei  f- interact  ion  *.**;*'  is  reduced  dramatically 
within  a  solid.  Aft  a  consequence,  nuclear  sei f- interact ion 
provides  the  dominant,  non-number  conserving  form  of  interaction 
at  one  or  a  set  of  lattice  sites.  A  further  consequence  of  the 
reduction  of  involves  an  uncoupling  of  the  motion  of  the 

center  of  mass  (r,r.)  from  variations  in  the  proton-neutron 
separation  (rn(F)  in  each  single  particle  D*  wave  function  v*d,ul 
through  the  separability  condition, 

f)  * 

Here,  ^#Uf(rsa)  varies  on  the  length  scale  of  the  electrostatic 
interaction,  while  the  spatial  variation  of  ifrnuc(rn.p)  occurs  over 
the  range  of  the  strong  interaction.  The  initial  state  D+  BBC 
component  becomes  occupied  when  sufficient  numbers  of  D+  ions 
become  indistinguishable  on  the  length  scale  of  the  electrostatic 
interaction  and  the  associated  single  particle  wave  functions 
tf.Ue(reo)  approach  the  single  particle  ttloeh  states  of  Eq.  2.  Eq. 
5,  then,  may  be  used  to  describe  the  many-body  electrostatic  wave 
function  for  the  D*  BBC  initial  state,  with  the  understanding  that 
each  Wannier  state  *t(r.)  is  multiplied  by  a  site  independent 
nuclear  wave  function  *nue,  which  depends  only  on  the  value  of 
rn.p  and  not  on  the  location  within  the  lattice  of  the  deuteron. 
Bounds  on  the  probability  of  nuclear  self-interaction  then  can  be 
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obtained  from  the  absolute  square  of  the  integral  of  the  multiply 
occupied  initial  state  product,  4>B(r1)  .  ...4s(rn),  multiplied  by  a 
more  compactly  distributed  product  of  final  state  wave  functions 
in  which  all  of  the  final  electrostatic  boson  functions  fall 
within  the  nuclear  force  range  rnuc.  From  the  "bosons  in  and 
bosons  out"  selection  rule  of  ref.  1,  we  are  led  to  explore  a 
reaction  in  which  four  initial  state  Wannier  functions  coalesce 
over  a  volume  Vnuc  =  4/3  ir  rnuc3,  forming  an  eight  fermion  "nuclear 
bag"  state  analogous  to  a  free  space  excited  8Be*  nucleus. 
Symbolically  we  will  refer  to  this  coalescence  reaction  as 

4(4,kp,r)  «*  6  (4,kp/rnuc)Bloch  ,  (6) 

where  0(4,kp,rnue)  represents  the  "nuclear  bag-like" 
configuration.  This  reaction  is  the  first  step  in  the  fusion 
reaction 

4D+Bloch  -  8Be‘Bl0Ch  -  2  AHe  (7) 
We  calculate  the  reaction  rate  from  the  Fermi  Golden  Rule7, 


ft  =  2n/h  2  <i  I  V I  f >  <f  I  V|  i>  6  (Ej-Ef)  ,  (8) 

where  |  i>  denotes  this  inital  D+  BBC  state.  The  sum  is  over 
final  states  |f>,  whe^e 


Hb  nb/2  Nl 

|f>  =(1/NbJ)1/2  2  ||  (1/NJ  s  exp ( ikp  R,)  • 

{ rm>  m=5  s=l 

^d(rx)^d(r2)^d(r3)  ^d(rj 


(9) 


where  <^d  is  a  compact  "defect"  function,  centered  at  the  single 
site  rfusion  where  fusion  takes  place  and  has  spatial  variation  in 
its  center  of  mass  coordinate  only  over  the  range  rnuc.  Also  ft  = 
h/  ( 2 »r )  ,  where  h  is  the  Planck  constant,  and  5  (Ej-Ef)  is  the  6- 
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function  in  which  Et  is  the  initial  total  energy  and  £,  is  the 
final  energy.  The  6-function  guarantees  conservation  of  eneigj 
between  initial  and  final  states.  The  operator  refers  to  the 
nuclear  self-interaction  potential. 

To  establish  a  lower  bound  on  the  required  BBC  concentration 
associated  with  the  observed  heating',  we  treat  tne  square  veil 
self-interaction  nuclear  potential  case,  in  which  ar.  isolated 
square  well  is  defined  by  an  upward  shift  in  kinetic  energy  ?.  cf 
each  D*  ion  over  the  characteristic  nuclear  dimension  r,ivi<  near 
^fusion *  The  volume  of  the  well  is  derived  frost  r,_.  -  ^.1  *  it 
cm,  appropriate  for  a  mass  8  nucleus.  This  is  equivalent  ti 
immediate  dissipation  of  a-particle  product  energy  *n  the 
reversed  reaction  within  the  region  where  interaction  takes 
place.  Consistent  with  these  approx iaat ions  we  replace  the  eight 
body  nucleon  problem  associated  with  the  four  defect  functions 
with  a  single  body  problem  in  which  the  four  initial  and  final 
state  electrostatic  functions  are  replaced  by  a  single  initial 
and  final  state  wdve  packet,  whose  characteristic  ewe  is 
determined  by  the  characteristic  length  scale  of  the  center  cf 
mass  motion  of  initial  and  final  states.  This  lc  equivalent  tc 
transfering  all  nuclear  energy  directly  into  alpha  particle 
motion  and  treating  all  nucloons  with  an  independent  nuclecr, 
model.  Also,  consistent  with  this  choice  of  model,  the  final 
density  of  states  p{ E)  *  dn/dE  is 

M°>  -  1/  (hu^J 

where  wnuc  =  (3h/2)  /(m,!^*),  »,  is  the  deuteron  mass,  and  t  -  -0 
is  the  energy  level  of  the  BBC  deuteron  band  in  the  lattice.  We 
choose  vanishing  values  for  all  wave-vectors  in  Eq.  5  since  the 
case  of  an  isolated  perturbation  applies  to  low  temperature. 
For  initial  and  final  states,  we  use  wave  packet  forms 

*p.ck.t.i  =  ( 2/ ( »aJ) ) 3/4  exp(-ry a2) 
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*  -  (V(2-Rw/ns  *  exp  (-31^/41^..^ 

where  r  is  measured  from  the  center  of  boson  Bass  in  each  unit 
cell  and  *  (V4»  Vfc...  )'  :.  The  overlap  integral  O(J)  is 

0(J)  -  (l6/3):’  (Vfc../V,.,);  * 


and 


ft  -  (16/3)’'  2*Vj/h  2m»R,„V(lh;)  c‘  .  (10) 

In  the  calculation  0(J)  receives  contributions  from  each  unit 
cell  in  the  crystal.  Using  the  heat  release  associated  with 
reaction  (7)  for  t.ie  magnitude  of  the  perturbation  ve  =  47.6  MeV, 
we  obtain  a  lower  bound  for  7  ■  2.8  x  10  ’  for  a  volumetric  power 
density  output  ftcyv.i l#  -  10  W/cm\  as  observed  in  ref  (2). 
is  the  unit  cell  volume  of  Pd.) 

Electrostatic  self-interaction  (the  coulomb  barrier) 

Inclusion  of  the  D*-D*  electrostatic  interaction  affects  the 
BBC  by  introducing  a  many-body  interaction  that  could  affect  the 
band  picture  from  which  the  BBC  is  derived.  The  band  picture 
remains  valid  provided  that  each  D*  (within  or  outside  the  BBC) 
remains  well -screened.  However,  a  breakdown  of  the  application 
of  the  band  picture  for  describing  nuclear  behavior  occurs  when 
the  timescale  r*1*'  associated  with  electrostatic  self-interaction 
for  multiply  occupied  Wannier  states  becomes  comparable  to  the 
timescale  associated  with  nuclear  self-interaction  rnuc.  The 
difference  between  r*1,c  and  rnuc  (=  the  inverse  of  the  nuclear 
reaction  rate)  determines  whether  or  not  the  electrostatic  and 
nuclear  interactions  are  separable,  and  thereby  determines  the 
appropriateness  of  the  form  assumed  for  each  mutually-non¬ 
interacting  boson  wave  function  in  the  nuclear  reaction 
calculation,  described  in  the  last  section. 
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EaCii  n-fold  occupation  component  of  the  BBC  is  a  single 
entity  whose  particle  character  is  evident  only  through  the 
existence  of  virtual  states  as  described  by  the  Wannier 
functions.  The  virtual  states  describe  a  transient  occupation  of 
a  unit  cell  by  an  integer  number  of  deuterons.  The  lifetime  of 
the  n-fold  multiple  occupation  virtual  state  is  determined  by  the 
self-interaction  energy  AEn,  except  for  the  single  occupation 
state.  Since  the  BBC  is  in  thermodynamic  equilibrium  "ith  non- 
BBC  deuterium  in  the  lattice,  aEx  must  be  small,  otherwise  BBC 
states  would  not  be  occupied.  AE2  is  the  difference  in  self 
interactio.i  between  single  and  double  occupation  of  a  unit  cell, 
and  can  be  evaluated  using  the  Wannier  state  wave  functions 
defined  by  Eq.  3: 

AE2  =  3/2  e2  J  J  d3r  d3r'  .  (11) 

I r-r ' j 

Using  a  =  0.5  Bohr,  we  find  aE2  =  81.6  eV.  Similarly  AE*  =  408 

eV.  The  corresponding  Planck  lifetimes  relocn  =  h/AEn  are  relec2  = 
5.1  x  10"17  s  and  tb1*c4  =  1.0  x  10‘1;  s.  These  electrostatically 
determined  lifetimes  compare  with  the  Planck  time  Tnuc*  =  h/vo 
10'22  s,  associated  with  self-induced  nuclear  interaction  that 
results  from  quadruple  occupation  (as  modeled  by  four 
distinguishable  single  particle  deuteron  wave  functions  which 
coalesce  within  a  common  nuclear  volume) .  The  large  difference 
in  electrostatic  and  nuclear  lifetimes  justifies  the  Born- 
Oppenheimer  separable  form  for  Vd.ut  used  in  the  last  section, 
which  leads  to  an  uncoupling  of  the  nuclear  and  electrostatic 
degrees  of  freedom.  The  reduction  in  charge  density  caused  by  D+ 
harmonic  motion  inside  the  unit  cell  thus  reduces  AE.lectrostatir  by 
five  orders  of  magnitude  from  its  free  space  value,  which  is 
comparable  to  AEnuc .  If  it  were  not  for  the  dominance  of  lattice 
interaction,  the  free  space  interaction  energy  would  apply  Born- 
Oppenheimer  separability  would  not  apply,  and  fusion  would  be 
prevented. 
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D*  BBC  Interactions  in  PdD 

The  last  section  describes  a  D+  BBC  volumetric  nuclear  self¬ 
interaction  that  is  calculable  using  the  Bloch  eigenstates  of  the 
system  and  that  approximately  conserves  crystal  momentum.  The 
reactions  do  not  refer  to  any  particular  value  of  x.  However, 
the  most  favorable  condition  for  D+  BBC  formation  is  x-l.  Under 
these  conditions  interaction  of  the  BBC  with  non-BBC  D+  may  be 
possible.  The  applicable  reactions  would  be 

Bloch  **  6ki  Bloch  (12a) 

6Li*virtual  +  2D+  -  ®Be*  -  2  alpha  .  (12b) 

The  maximum  rate  for  this  reaction  chain  would  have  the  same  form 
as  that  of  Eq.  7,  but  with  a  c3  concentration  dependence  instead 
of  a  c‘  dependence,  namely 

n  =  ( 16/3 ) 3/2  2 jtV02/!^  2mBRnuc2/(3*2)  Vnuc/Vhar  c3  .  (13) 

An  alternate  mode  of  energy  release  might  be  by  successive 
nuclear  interaction  scattering  reactions  with  non-BBC  D+.  These 


reactions  could 

be  of  the  form 

2D  BIoch  ** 

‘He*Bloch 

(14a) 

* 

virtual 

D  -  ‘He**  +  Dscat 

(14b) 

‘He**  +  D  -» 

‘He*"  +  Dacat 

(14c) 

‘He**'  *  +  D 

-  ‘He  +  Dseat 

,  ( 14d) 

where  4He**  and  4He**"‘  are  nuclear  bag  states  of  lower  total 
energy  than  the  D+  BBC  and  the  Dacat  are  recoil  scattered  non-BBC 
deuterons.  Once  reaction  (14b)  occurs,  reaction  (14a)  becomes 
energetically  blocked.  Reactions  (14)  may  be  the  most  favorable 


mode  of  decay.  Coupling  between  the  chemically  bonded  D  and  the 
lattice  could  result  in  momentum  transfer  to  the  lattice  as  a 
whole.  The  result  could  be  heat  release  without  any  high  energy 
particle  generation. 

(a)  Permanent  address:  Bendix  Field  Engineering  Corp.  ,  Oxon  Hill, 
MD  20745 
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